Lateral Casimir force beyond the Proximity Force Approximation 
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We argue that the appropriate variable to study a non trivial geometry dependence of the Casimir 
force is the lateral component of the Casimir force, which we evaluate between two corrugated metal- 
lic plates outside the validity of the Proximity Force Approximation (PFA). The metallic plates are 
described by the plasma model, with arbitrary values for the plasma wavelength, the plate sepa- 
ration and the corrugation period, the corrugation amplitude remaining the smallest length scale. 
Our analysis shows that in realistic experimental situations the Proximity Force Approximation 
overestimates the force by up to 30%. 



Considerable experimental progress has been 
achieved in the measurement of the Casimir 

force, opening the way for various applications in 
nano-science particularly in the development of 
nano- or micro-electromechanical systems (NEMS or 
MEMS). Calculations are much simpler in the original 
Casimir geometry of two plane plates which obeys 
a symmetry with respect to lateral translations and 
thus allows to derive the expression of the Casimir force 
from the reflection amplitudes which describe specular 
scattering on the plates 4| . 

More general geometries open a far richer physics 
with a variety of extremely interesting theoretical pre- 
dictions ^] . Up to now the experimental studies of the 
effect of geometry have been restricted to simple configu- 
rations which can be calculated with the help of the Prox- 
imity Force Approximation (PFA). This approximation 
is essentially equivalent to an averaging over plane-plane 
geometries and its result can be deduced from the force 
known in this geometry 0. For example, it allows to 
evaluate the force between a plane and a sphere pro- 
vided the radius R of the sphere is much larger than the 
mirror separation L. ft is also valid for the descrip- 
tion of the effect of roughness when the wavelengths as- 
sociated with the plate deformations are large enough @ . 
However PFA relies heavily on assuming some additivity 
of Casimir forces which is known to be generally not valid 
except for very smooth geometrical perturbations ;9j . 

The aim of the present paper is to study a configu- 
ration allowing a new test of QED theoretical predic- 
tions outside the PFA domain and independent of those 
already performed in the plane-plane geometry. The 
idea is to look for the lateral component of the Casimir 
force which appears, besides the usual normal compo- 
nent, when periodic corrugations with the same period 
are imprinted on the two metallic plates. This configu- 
ration contrasts with other ones, for example the normal 
Casimir force in the plane-sphere geometry or roughness 
corrections to it. There PFA can also be invalid, but 
this leads only to small corrections of the dominant nor- 
mal Casimir force, which do not seem accessible experi- 
mentally at the moment. The lateral component of the 



Casimir force has recently been measured and analyzed 
within the PFA We find for experimentally re- 

alizable parameters that PFA overestimates the force by 
as much as 30%, which should allow for a rapid experi- 
mental check of its validity. 

The lateral Casimir force between corrugated surfaces 
has been analyzed outside the PFA domain for perfect 
reflectors 12] where interesting results were obtained for 
arbitrary values of the ratio Xc/L of the corrugation 
wavelength Ac to the interplate distance L. However, the 
assumption of perfect reflections can only be valid in the 
limit of large distances where the lateral force tends to 
become too weak to be measurable. Experimental condi- 
tions allowing the lateral force to be measured correspond 
to separations of a few hundred nanometers ||lOj, that 
is of the same order of magnitude as the plasma wave- 
length Ap associated with the metallic plates, so that 
they cannot be treated even as approximately perfect re- 
flectors In contrast, the calculations presented here 
allow to predict the lateral force in this domain. 

In this letter, we calculate the lateral force for metallic 
plates modeled by the plasma model with arbitrary val- 
ues of L, Ac and Ap. We use the perturbative approach 
that we developed for analyzing the effect of stochastic 
roughness on the normal Casimir force [11 m. This 
technique is valid as long as the corrugation amplitude 
a remains the smallest of length scales a <^ L,Ac,Ap. 
As this condition does not depend on the relative mag- 
nitudes of the three other length scales, it allows us to 
study situations beyond the PFA where the lateral force 
is experimentally accessible. The result will be expressed 
in terms of a non linear susceptibility function, calcu- 
lated from non-specular scattering amplitudes flq asso- 
ciated with corrugated metallic surfaces. This function, 
obtained within the scattering approach [T^, can itself 
be considered as a new QED theoretical prediction to be 
compared with forthcoming experiments. 

The surface proflles of the two parallel corrugated 
plates are defined by two functions /ii(r) and /i2(r), with 
r = {x,y) the lateral position along the surfaces of the 
plates. Both distributions hi and /i2 have zero spatial 
averages (hj) = 0,j = 1,2, and they are counted as 
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positive when they correspond to local length decreases 
below the mean value L. The corrugated surfaces are 
assumed to be static, so that the field frequency w is pre- 
served by scattering. In contrast, the lateral wavevector 
components k = {kx,ky) as well as the polarization of 
the field are modified. Scattering on the plate j — 1,2 
is thus described by non-specular reflection amplitudes 
TZj-pp' (k, k', w), where k, k' represent the lateral wavevec- 
tors of the input and output fields and p and p' their po- 
larizations, TE for transverse electric and TM for trans- 
verse magnetic (more detailed definitions in [T^lisl l'). 

For the purpose of the present paper, the non-specular 
refiection amplitudes have to be developed up to the first 
order in the deviations hj from fiatness of the two plates 
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(1) 



The first line in this equation represents the zero-th or- 
der term with respect to corrugation, that is also the 
specular reflection on a flat plate (wit h r., :n(k, u) the or- 
dinary specular reflection amplitude |l7j). The second 
line describes the first-order correction proportional to 
the Fourier component Hj(k — k') of the profiles hj{r), 
this Fourier component being able to induce a modifica- 
tion of the field wavevector from k to k'. 

We then compute the correction of the Casimir energy 
SEpp induced by the corrugations. At the lowest order, 
we have to use the scattering approach at second or- 
der in the corrugations, keeping only the crossed terms 
of the form Hi H2 which have the ability to induce lat- 
eral forces. This means that the sensitivity function ob- 
tained below depends on the crossed correlation between 
the profiles of the two plates, in contrast to the func- 
tion which was calculated for describing roughness cor- 
rection in 0, . The latter were depending on terms 
quadratic in Hi or H2 , and their evaluation required that 
second order non specular scattering be properly taken 
into account. Here, first order non specular amplitudes 
evaluated on both plates are sufficient. 

Assuming for simplicity that the two plates are made 
of the same metallic medium, so that explicit reference to 
the index j may be omitted in the reflection amplitudes 
from now on (otherwise the result is given by a trivial 
extension), we obtain the second-order correction 
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The non linear response function Qik) is given by 

e(k) = -n/ ^/(^^.^.-.(O (3) 
g-(K'+«)L^^,^(k/^ k; i)Rpp, (k, k'; 



6k', k = 



P ,P 



(l-rp'(k')2e-2K'L)(i_^^(k) 



2 „-2kL 



Note that the integral over real frequencies lo has been re- 
placed by an integral over imaginary frequencies ^, while 



the round-trip propagation factor between the plates is 
now expressed as exp(— 2/tL), with k = a/Zc^ _j_ ^ For 
isotropic media, symmetry requires the response function 
C7(k) to depend only on the modulus fc = |k|. 

Experiments with corrugated plates ^IQj were corre- 
sponding to the simple case where uniaxial sinusoidal 
corrugations are imprinted on the two plates along the 
same direction, say the y direction, and with the same 
wavevector k = 27r/Ac 



h\ — a\ cos {kx) 



h2 — a2 cos {k{x + b)) (4) 



The energy correction thus depends on the lateral mis- 
match h between the corrugations of the two plates, which 
is the cause for the lateral force to arise. Replacing the 
ill-defined {2iTf 5^'^\Q) by the area A of the plates, we 
derive from (jSJ 



6Epp = A'^cos{kb)g{k) 
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The result of the PFA can be recovered from Eq. Q 
as the limiting case fc ^ 0. PFA indeed correspond to 
long corrugation wavelengths, that is also nearly plane 
surfaces so that the Casimir energy is obtained from the 
plane-plane energy Epp{z) by averaging the 'local' dis- 
tance z — L — hi ~ h2 over the surface of the plates |^. 
Expanding at second order in the corrugations, we thus 
obtain the following dependence of the energy (as before 
we disregard terms in and because they do not 
depend on the relative lateral position) 

SEpp = ^cos{kb)^ , [PFA] (6) 
2 dz-^ 

This reproduces lO for large values of Ac because the 
response function g{k) satisfies the condition 



lim Ag{k) = '^-^ 
k^o dz^ 
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This property is ensured, for any specific model of the 
material medium, by the fact that Q is given for k —^ 
by the specular limit of non specular refiection ampli- 
tudes (more details in ^3)- For arbitrary values of k, 
the deviation from PFA is described by the ratio 
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In the following, we discuss explicit expressions of this 
ratio p obtained from the general result Q in the par- 
ticular case of metallic mirrors described by the plasma 
model. The dielectric function is thus given by e = 
1 -|- cjp/^^^ with the plasma wavelength and frequency 
related by Ap = 2^0/ up. For the numerical examples, 
we will take Ap = 136nm, corresponding to gold covered 
plates. The non-specular coefficients appearing in Eq. Q 
are then evaluated from Ref. ISi]. The resulting function 
p is plotted on Fig. ^ as a function of k, for different 
values of the separation distance L. 
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FIG. 1: Variation of p versus k with Ap = 136nm and for L = 
50nm (dotted line), lOOnm (dash-dotted line), 200nm (solid 
line) and 400nm (dashed line). 
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FIG. 2: Lateral force amplitude for the plane-sphere setup, as 
a function of k, with figures taken from ^XJ\. The experimental 
value k — 0.0052nm~^ is indicated by the vertical dashed line. 

As expected, the PFA (p = 1) is a better approxima- 
tion at short distances. For L =50nm for example, the 
approximation is correct in the range k < O.Olnm"^ (i.e. 
Ac > 628nm) covered by the plot in Fig. ^ For larger 
distances, PFA is found to overestimate the energy cor- 
rection and, therefore, the lateral force. For parameters 
L =200nm and Ac = 1.2/zm {k = 0.0052nm"i) close to 
the experimental figures of Ref. ^3 1 we find the lateral 
force smaller by a factor p — 0.84. In other words, the 
correction expected with respect to PFA calculations is 
of the order of 16%. 

For still larger values of fc, the functions G{k) and p{k) 
decay exponentially to zero. When the momentum trans- 
fer k is much larger than 1/L, the function 6k',k'-k in 
Eq. 10) is proportional to the exponentially small factor 
exp(— /cL). If we also assume that fcAp ^ 1, we find 
G{k) ~ ak exp(— /cL) where the parameter a now de- 
pends on Ap and L only. This is in striking contrast 
with the behavior of the response function for stochastic 
roughness, which grows linearly with k for large k due to 
the contribution of the second-order reflection coefficients 
Il5l. These coefficients do not contribute to the second- 



order lateral effect, which is related to two first-order 
non-specular reflections at different plates, separated by 
a one-way propagation with a modified momentum of 
the order of k. The resulting propagation factor is, in the 
large-/c limit, exp(— kL) « cxp(— fcL), thus explaining the 
exponential behavior. 

For intermediate values of fc, Q{k) is also smaller than 
the roughness response function, although no simple an- 
alytical expression is available in the general case. Thus, 
the exact calculation decreases the PFA value for the lat- 
eral effect, in contrast to the case of stochastic roughness 
correction where they increase the PFA value [l4L Il5| . 

We now compare the theoretical expression of the lat- 
eral Casimir force to realistic experiments and therefore 
consider the plane-sphere (PS) geometry 1^ rather than 
the plane-plane (PP) geometry as before. The PS result 
can be derived from the PP one by using PFA with va- 
lidity conditions more easily met than those which would 
be necessary for computing the corrugation effect. As a 
matter of fact, PFA may account for the curvature of the 
sphere as soon as the radius R of the sphere is larger than 
the distance R ^ L. Then any interplay between curva- 
ture and corrugation is avoided provided that RL ^ A^ . 
These two conditions are met in the experiment reported 
in 0, where R = lOO/xm, Ac = 1.2/im and L ^ 200nm. 
PFA is likely a very good approximation as far as cur- 
vature is concerned whereas, for the corrugation effect 
itself, a deviation from the PFA is expected since L is 
not so much smaller than Ac- 

Applying the PFA to describe the PS configuration, 
we obtain the energy correction 6Eps between the sphere 
and a plane at a distance of closest approach L as an in- 
tegral of the energy correction SEpp in the PP geometry 

dEpsiL, b) = '^-^^^5Epp{L\ b) (9) 

Then the lateral force is deduced by varying the energy 
correction with respect to the lateral mismatch b between 
the two corrugations. This gives the lateral Casimir force 
in the PS geometry as an average of 

Fn{L^b)^-^^Ep,{LM^ jy-^^F^^^{L',b) 
Flf{L,b) = -^Epp{L,b) (10) 
leading eventually with Eq. ISJ to 

F]ifs' = 7raia2fci?sin(fc6) / dL'g{k,L') (11) 

J oo 

The force attains a maximal amplitude for sin(fc6) — ±1, 
which is easily evaluated in the PFA regime fc ^ where 
Q{k) does not depend on /c, so that i^pg' scales as k. How- 
ever, as k increases, the amplitude increases at a slower 
rate and then starts to decrease due to the exponential 
decay of Q{k). Thus, for a given value of the separation 
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FIG. 3: Plane-sphere geometry: _F]!)'s/(aia2 sin(fc6)) as a func- 
tion of L, with k = 0.0052nm~^. The solid line is the exact 
result obtained in this paper. Other lines represent various ap- 
proximations: PFA (dotted-dashed) and PFA combined with 
the perfectly-reflecting (dashed) and plasmon (dotted) limits. 

L, the lateral force reaches an optimum for a corrugation 
wavelength such that kL is of order of unity, which gen- 
eralizes the result obtained for perfect reflectors in ,12] . 
In Fig. 13 we plot the force Fpg* (for sin(fc5) = 1) as 
a function of fc, with figures taken from the experiment 
of Ref. We also use the values ai = 59nm and 

02 = 8nm of the amplitudes for measuring the force as 
in [lOf , reminding however that our calculations are valid 
in the perturbative limit ai, a2 0. 

The plot clearly shows the linear growth for small k 
as well as the exponential decay for large k. The maxi- 
mum force is at /c = 0.009nm~^ so that kL ~ 2. The 
experimental value k — 0.0052nm^^ is indicated by the 
dashed line in Fig. 12 and the force obtained as 0.20pN, 
well below the PFA result, indicated by the straight line 
and corresponding to a force of 0.28pN. This last value 



is very close to the number calculated by Ref. [lOj, if 
we discount the correction due to higher order terms be- 
yond the second order approximation considered in the 
present paper. Precisely, Ref. 10] finds a force of 0.32pN 
ai L — 221nm, with a relative correction due to higher 
powers of 1.21. Discounting this factor, the second-order 
force should be 0.26pN, which overestimates the exact 
result by a factor of the order of 30%. 

In order to discuss the dependence of the force versus 
the distance L, we plot J^pg*/(aia2 sin(fc6)) with log-log 
scales on Fig. 13 with Ac = 1.2//m (solid fine). For very 
small values of L, we find the power law expected for 
the plasmon limit L ^ Ap "T^ combined with PFA (dot- 
ted line) . We also show the PFA result for arbitrary \p/L 
(dotted-dashed line) and the l/L^ dependence resulting 
from the combination of PFA and perfectly-reflecting as- 
sumption (dashed line). Fitting in the range between 
150 and 300 nm, we find the power law L^^-^, in agree- 
ment with the experimental result fTol|. But Fig. 13 shows 
that the force decays faster for larger values of L, leading 
ultimately to an exponential decay. 

These results indicate the way to be followed to per- 
form an accurate comparison between theory and experi- 
ment in a configuration where geometry plays a non triv- 
ial role, i.e. beyond the PFA. Such a comparison should 
be pursued by looking at the lateral force between cor- 
rugated metallic plates. This problem can neither be 
addresse d by calculations performed with the help of 
PFA [13 ' '^O'^ with the assumption of perfect refiec- 
tion jl4|. But it can be tackled by the technique pre- 
sented in the present letter which is valid for arbitrary 
relative values of L, Ac and Ap, provided that the corru- 
gation amplitude is smaller than the other length scales. 
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